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摘 要 : 利用 多 元 线性 回归 模型 (MLR) 、 主 成 分 逐步 回归 分 析 模 型 (PCSR ) 、 克 里 金 插值 法 (Kriging) ,将 青海 湖 流 
域 及 周边 区 域 0.25° 分 辨 率 的 TRMM 3B43 降水 数据 降 尺 度 至 0.01° 分 辨 率 , 并 选取 研究 区 范围 内 20 个 气象 站 点 的 实 
测 降 水 数据 ,利用 相关 系数 (CC) \ 均 方 根 误差 (RMSE) 和 相对 偏差 (Bias) 对 降 尺 度 结果 进行 评价 ,以 此 选取 更 适用 
于 研究 区 的 降 尺度 方法 。 结 果 表 明 :(1) 基于 TRMM 及 3 种 降 尺 度 方法 获取 的 研究 区 降水 空间 分 布 具 有 一 致 性 ,年 
均 及 春 、 夏 .秋季 的 降水 量 均 表 现 为 北部 高 ,西部 及 西北 部 低 , 而 冬季 降水 量 表现 为 南部 与 西北 部 高 ,中 部 低 -(2) WE 
究 区 降水 量 随 着 海拔 的 增高 ,以 3800m 为 界 整体 上 呈现 先 升 高 后 降低 的 趋势 。(3) 精度 评价 的 结果 表明 ,年 尺度 上 
的 Kriging 精 度 表 现 最 好 ;在 空间 上 TRMM 及 3 种 降 尺 度数 据 在 东部 地 区 精度 最 优 。 在 季 尺 度 上 ,数据 精度 表现 为 
PCSR>Kriging>TRMM>MLR ;在 月 尺度 上 ,PCSR 数据 精度 最 优 .(4) 海拔 对 研究 区 内 TRMM 及 3 种 降 尺 度数 据 的 影 


响 较 小 ,但 随 着 海拔 的 升 高 ,遥感 数据 逐渐 出 现 低估 降水 的 现象 ,其 可 能 原 


因 在 于 降水 的 低估 与 微波 降水 率 反 演 时 


对 对 流 性 降水 的 低估 有 关 。 综 合 降水 空间 分 布 一 致 性 分 析 与 精度 评价 ,认为 PCSR 最 适合 于 青海 湖 流域 及 周边 区 


域 的 TRMM 3B43 降水 数据 降 尺 度 方法 。 
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降水 是 气候 系统 中 的 一 个 重要 参数 ,是 水 循环 
的 重要 组 成 部 分 ,准确 的 降水 数据 对 理解 区 域 乃 至 
全 球 气 候 变 化 和 水 循环 过 程 具 有 重要 意义 "*。 H 
前 ,降水 数据 的 获取 方式 包括 地 面 气象 站 观测 Ba 
古 雷 达 和 卫星 遥感 探测 ”“' 。 在 这 3 种 方法 中 ,地 面 
观测 获取 的 降水 量 最 为 准确 ,高 密度 分 布 的 雨量 站 
理论 上 可 以 获取 准确 的 降水 空间 分 布 特征 ,然而 在 
高 海拔 山区 流域 ,由 于 观测 站 点 稀少 且 降 水 空间 差 
异性 大 ,准确 的 获取 这 些 区 域 降水 空间 分 布 信息 至 
今 是 一 个 挑战 “。 降 十 雷达 通过 回 波 强度 对 地 面 
进行 降水 预测 ,可 以 减 小 气象 站 点 以 点 代 面 的 误 
Fe ,但 其 约束 条 件 较 多 TE PERE OP), TU 
感 探 测 具有 和 履 盖 范围 广 、 时 空 分 辩 率 高 和 约束 条 件 
少 的 特点 ,能够 弥补 地 面 观测 的 不 足 ”'"。 随 着 遥 
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感 技术 的 发 展 ,不 断 出 现 的 卫星 降水 产品 在 气象 水 
文人 研究 领域 得 到 了 广泛 应 用 ,如 CMORPH"、 
GPCP' TRMM™"" PERSIANN-CDR "I 以 及 GPM 
IMERG!” 等 。 

然而 , 现 有 的 遥感 降水 产品 普遍 分 辩 率 低 , 难 
以 满足 山区 小 流域 的 研究 要 求 , 因 此 对 各 感 降水 产 
品 的 降 尺 度 研 究 逐 渐 成 为 热点 。Elnashar 等 ”利用 
Google Earth Engine 和 Google Cloud Computing 对 
TRMM 月 降水 量 进行 降 尺 度 分 析 ; 范 雪 苓 等 ”基于 
TRMM 数据 在 天 山 山区 构建 主 成 分 -逐步 回归 模 
型 ,结果 表明 降 尺 度 后 的 数据 精度 显著 优 于 原始 
TRMM 3B43 数据 ; 张 寒 博 等 3 对 GPM IMERG 和 
TRMM 3B43 数 据 进行 空间 隆 尺度 ,以 此 对 华中 地 区 
进行 降水 时 空 变 化 分 析 和 干旱 监测 ; 李 豪 等 所 采用 
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MGWR 与 Kriging 相 结合 的 方法 ,建立 了 MGCWRK 模 
型 进行 降 尺 度 ,结果 表明 降 尺 度 处 理 不 能 提高 研究 
区 数据 精度 ; 徐 彬 仁 等 ”基于 随机 森林 算法 对 青藏 
高 原 TRMM 降水 数据 进行 了 空间 统计 降 尺 度 研 究 ; 
Ghorbanpour 等 "站 利用 5 种 降 尺 度 方法 对 马尔 米 亚 
湖 流域 进行 降 尺 度 评 估 。 类 似 的 降 尺 度 研究 还 有 
很 多 ,但 总 体 上 并 没有 一 个 普 适 的 降 尺 度 标准 。 

青海 湖 流 域 位 于 青藏 高 原 东 北部 ,是 重要 的 生 
态 屏 障 “”。 流 域内 国家 气象 局 公开 数据 的 气象 站 
点 仅 有 1 个 ,对 青海 湖 流域 的 水 文 研究 造成 很 大 限 
制 ?”。 张 乐 乐 等 "在 青海 湖 流 域 及 周边 地 区 通过 
精度 验证 评价 了 TRMM 3B42V7, CMORPH , PERSI- 
ANN, PERSIANN- CDR 遥感 降水 资料 ,结果 表明 
TRMM 数 据 记 录 的 降水 量 最 为 准确 ,精度 最 好 ,更 适 
用 于 青海 湖 流域 及 周边 地 区 。 然 而 ,TRMM 数据 由 
于 空间 分 辩 率 低 (0.25°x0.25°) ,难以 满足 小 流域 尺 
度 水 文 过 程 研 究 需 求 。 因 此 ,有 必要 通过 降 尺 度 获 
取 更 高 分 辨 率 的 降水 数据 。 

因 青 海 湖 流域 气象 站 点 稀少 ,本 人 研究 将 研究 范 
围 扩大 到 包含 周围 20 个 气象 站 点 的 区 域 ,并 以 此 作 
为 研究 区 ,选用 多 元 线性 回归 模型 (Multivariable 
Linear Regression, MLR) 、 主 成 分 -逐步 回归 模型 
(Principal Component Stepwise Regression , PCSR ) 和 
克 里 金 插值 法 (Kriging) 对 TRMM 3B43 数据 进行 降 
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尺度 ,并 利用 研究 区 内 20 个 气象 站 点 的 实测 数据 ， 
从 年 、 季 、 月 尺度 对 降 尺 度 结 果 进 行 精度 评价 。 本 
人 研究 的 主要 目的 是 选 出 最 合适 的 青海 湖 流域 
TRMM 3B43 数据 降 尺 度 方法 ,为 流域 开展 高 分 辩 率 
生态 .水文 过 程 模拟 提供 数据 支撑 。 


1 研究 区 与 数据 方法 


1.1 研究 区 概况 

青海 湖 流域 位 于 36"15'~38?20'N ,97°50'~101°20'E 
之 间 ,流域 面积 约 29661 ko ,湖泊 面积 约 为 4400 km’, 
青海 湖 流域 南 至 青海 南山 , 西 至 阿木 尼 尼 库 山 , 东 
至 日 月 山 , 北 至 大 通 山 ,形成 一 个 完整 的 封闭 型 高 
原 内 陆 盆 地 。 整 个 流域 呈 西 北 -东南 走向 ,地 势 西 
北 高 东南 低 光 2 。 因 青海 湖 流域 气象 站 点 稀少 , 研 
究 区 划 定 范围 扩大 到 包含 周围 20 个 气象 站 点 的 区 
域 ,经 纬度 范围 为 35°~39°N,97°~103°E( 图 1)。 研 
究 区 地 形 以 高 原 为 主 , 地 势 整体 西高 东 低 ,海拔 
1330~5663 m, ,西北 和 西南 整体 海拔 均 在 4000 m 以 
上 。 人 研究 区 降水 量 特征 表现 为 季节 分 布 不 均 , 夏 季 
降水 量 占 主导 地 位 且 较 为 集中 ,其 次 是 春秋 季 , 冬 
BR, 
1.2 数据 来 源 与 预 处 理 

TRMM 3B43 降水 数据 来 源 于 NASA 网 站 (https: 
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图 1 研究 区 位 置 及 气象 站 点 分 布 


Fig. 1 Location ofthe study area and distribution of meteorological sites 


202301.00153v1 


chinaXiv 


1708 干 时 


ChinaXiv 合 作 期 刊 


39 卷 


//search.earthdata.nasa.gov) , 空间 分 辨 率 为 0.25°x 
0.25° ,时 间 分 辩 率 为 月 尺度 ,覆盖 范围 在 全 球 南北 
纬 50" 之 间 ,时间 为 1979 年 至 今 , 本 研究 选取 的 数据 
年 份 为 2010—2019 年 。TRMM 数据 为 NetCDF 格 
式 , 处 理 时 首先 要 借助 ArcGIS 进行 格式 .投影 转 
换 。TRMM 3B43 原始 降水 数据 的 单位 为 mm'h-, 乘 
以 每 个 月 的 小 时 数 即 可 得 到 月 降水 量 。 

DEM 数据 采用 SRTMDEM 90 m 分 辨 率 原始 高 
程 数据 ,数据 来 源 于 地 理 空 间 数据 云 (http://www.gs- 
cloud.cn/search) 。 坡 度 和 坡 向 数据 由 DEM 数据 生 
成 。 并 通过 裁剪 、 重 采样 获得 研究 区 0.25°* 和 0.01° 
分 辩 率 的 高 程 数据 ,并 基于 ArcGIS 水 文 分 析 工具 得 


PCSR 方 法 可 以 将 多 个 影响 因子 转化 成 几 个 独 
立 因 子 进行 分 析 。 在 保证 信息 最 大 化 的 前 提 下 对 
原来 的 影响 因子 进行 处 理 , 并 在 处 理 过 程 中 剔除 不 
显著 的 解释 变量 ,以 此 来 提高 研究 效率 并 使 信息 更 
为 精确 2 。 本 文通 过 构建 PCSR 模型 建立 TRMM 数 
据 与 影响 因子 的 方程 ,进行 降 尺 度 空间 转换 ,模型 
构建 如 下 : 
Y=(b+b X, +b,X,+b Xs +:*")+Kk (2) 
式 中 :7Y 为 TRMM 降 水量 降 尺 度 后 的 数据 ;5b~b; 表 示 
主 成 分 回归 系数 ;~ 名 表示 主 成 分 回归 因子 ; (b+ 
b,X\+b.X,+b:X3+ +++ ) 表 示 由 PCSR 回归 得 到 的 降水 


到 了 研究 区 0.25° 和 0.01° 分 辨 率 的 坡度 和 坡 向 

NDVI 数 据 采 用 MODIS 数 据 MOD13A3 产 品 , 数 
据 来 源 于 NASA 网 站 (https://search.earthdata.nasa. 
gov) ,空间 分 辨 率 为 1 km, 时 间 分 辩 率 为 月 尺度 , 数 
据 下 载 年 份 为 2010 一 2019 年 。 处 理 时 首先 运用 
ENVI 进 行 图 像 拼 接 .异常 值 处 理 、 格 式 转换 ,再 在 
ArcGIS 中 进行 投影 转换 和 研究 区 裁剪 。 

气象 站 点 观测 降水 数据 采用 中 国 气象 数据 网 
(http://data.cma.cn) ,在 本 研究 区 内 共有 20 个 气象 观 
测 站 点 ,分 别 为 托 勒 .野牛 沟 . 祁 连 、 德 令 喻 、 刚 察 、 
门 源 .马兰 AB SEA ET 贵 德 . 民 和 、 兴 海 、 贵 
南 .同仁 ` 茶 卡 SHR .武威 ,永昌 、 山 丹 ( 图 1)。 本 
研究 选取 的 降水 观测 数据 时 间 尺 度 与 TRMM 数据 
一 致 , 观 测 数 据 的 时 间 分 辩 率 为 日 ,通过 佬 加 可 以 
获取 月 .季节 和 年 太 度 的 降水 数据 ,研究 期 内 数据 
无 缺 测 。 
1.3 研究 方法 

降 尺 度 方法 选用 多 元 线性 回归 模型 (MLR) 、 主 
成 分 -逐步 回归 模型 (PCSR) 和 普通 克 里 金 插值 法 
(Kriging)。 其 中 Kriging 方 法 可 以 在 ArcGIS 软件 中 
直接 实现 。MLR 通过 2 个 或 2 个 以 上 的 影响 因素 作 
为 自 变 量 来 解释 因 变 量 的 变化 ,以 此 建立 影响 因子 
与 自 变 量 的 回归 模型 。 降 水 的 形成 受到 多 种 因素 
的 影响 ,本 文选 取经 度 .纬度 ,高程 .坡度 . 坡 向 和 
归 一 化 植被 指数 (NDVI) 作 为 影响 因子 ,建立 与 
TRMM 降水 数据 的 回归 模型 ,一 般 方程 如 下 : 

y=ag tax +a, x, +a% +a% +a; x; taoexe tk (1) 

式 中 :y 为 TRMM BEJK ZE ; ao~as R AN EIJ AL ; 
xi~xe 表 示 回 归 因 子 ;k 表 示 降 水 残 差 值 。 


回归 值 ;K 表 示 降 水 残 差 值 。 

具体 步骤 如 下 :将 DEM 数据 .经 纬度 .NDVI 处 
理 到 相同 0.25° 分 辩 紊 下 ,运用 Matlab 软件 将 各 数据 
所 有 栅 格 中 心 像 元 点 的 数值 进行 汇总 ;通过 Matlab, 
SPSS 计 算得 到 回归 系数 ,汇总 得 到 0.25° 分 辨 率 降 
水 预测 回归 值 和 降水 残 差 ;采用 普通 克 里 金 法 对 降 
水 残 差 进行 插值 ,得 到 0.01" 分 辩 率 残 差 数据 ;将 
0.25? 分 辨 率 下 得 到 的 回归 系数 代入 到 0.01? 分 辩 率 
下 的 回归 因子 中 ,并 加 上 高 分 辩 率 残 差 数据 ,最终 
得 到 0.01° 分 辩 率 下 TRMM 降水 量 降 尺 度数 据 , 实 现 
降 尺 度 的 转换 。 
1.4 精度 评价 

选择 2010 一 2019 年 年 均 . 季 均 和 月 份 的 降 尺 度 
数据 ,通过 ArcGIS 提 取 研 究 区 范围 内 20 个 气象 站 点 
的 值 ,与 降水 实测 值 进行 精度 评价 ,采用 相关 系数 
(CC) 、 均 方 根 误差 (RMSE) 和 相对 偏差 (Bias) 评 价 
降 尺 度 结 果 ””Y。 其 中 CC 的 取 值 范围 为 [-1,1],CC 
值 越 接 近 1, 说 明 格 点 数据 与 观测 降水 数据 的 一 致 
性 越 高 。RMSE 的 取 值 范围 为 [0, +%), 值 越 接 近 0， 
表示 格 点 降水 数据 越 准 确 。Bias 的 取 值 范围 为 [-1， 
1] 值 , 取 值 越 接近 0 ,表示 格 点 降水 数据 与 观测 降水 
数据 偏离 程度 越 小 。 


2 结果 与 分 析 


2.1 青海 湖 流域 降水 空间 分 布 特征 

利用 MLR .PCSR 和 Kriging 方 法 对 TRMM 数据 
降 尺 度 ,并 通过 不 同时 间 尺 度 处 理 , 得 到 青海 湖 流 
域 2010 一 2019 年 年 平均 和 春 、 夏 、 秋 、 冬 4 个 季节 平 
均 降 水 量 的 空间 分 布 ( 图 2)。 图 2a~ 图 2d 给 出 了 
TRMM 原始 数据 以 及 3 种 降 尺 度 方法 得 到 的 青海 湖 
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图 2 青海 湖 流域 降 尺 度 结果 空间 分 布 
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Fig. 2 Spatial distribution map of downscaling results in Qinghai Lake Basin 
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流域 范围 内 年 平均 降水 量 空间 分 布 。 从 结果 来 看 ， 

经 过 降 尺 度 之 后 得 到 的 年 平均 降水 量 与 原始 TRMM 
数据 的 空间 分 布 一 致 ,降水 量 高 值 区 主要 分 布 于 青 
海 湖北 部 , 低 值 区 主要 集中 在 青海 湖西 部 及 西北 
部 。 从 降水 量 的 取 值 范围 来 看 ,Kriging 降 尺度 降水 
量 最 大 值 与 TRMM 数据 保持 一 致 ,MLR 降 尺 度 降 水 
量 最 大 值 略 高 于 TRMM 数据 ,而 PCSR 降 尺 度 降水 
量 最 大 值 略 低 于 TRMM 数据 ,此 外 ,3 种 降 尺 度数 据 
的 降水 量 最 小 值 均 低 于 TRMM 数据。 

由 图 2e~ 图 2t 可 知 , 经 过 降 尺 度 之 后 得 到 的 季 
均 降 水 量 与 TRMM 数据 的 空间 分 布 一 致 。 其 中 ， 
ER .秋季 空间 分 布 与 年 均 降 水 空间 分 布 类 似 , 而 
冬季 降水 量 高 值 区 分 布 在 青海 湖南 部 与 西北 部 , 低 
值 区 分 布 在 青海 湖 流域 中 部 。 

为 了 进一步 给 出 青海 湖 流 域 及 周边 区 域 不 同 


海拔 降水 量 的 差异 ,按照 500 m 间 隔 在 研究 区 划分 
出 了 9 个 海拔 梯度 ,并 统计 了 每 个 海拔 梯度 的 年 平 
均 降水 量 ( 图 3)。 其 中 ,5300~5663 mn 海拔 范围 内 空 
间 较 少 , 因 此 将 最 后 2 类 合并 为 4800~5663 m 一 
级 。 结 果 表 明 , 随 着 海拔 的 增高 ,研究 区 年 平均 降 
水 量 整体 上 呈现 先 升 高 后 降 休 的 趋势 。 在 1330~ 
2800 m, 降 水量 随 着 海拔 的 升 高 不 断 增 加 ,在 2800~ 
3300 m, 降水 量 随 着 海拔 的 升 高 有 小 幅 下 降 ,在 
3300~3800 m 海 拔 范围 内 ,降水 量 随 海 拔 的 升 高 而 
增加 并 达 最 大 值 , 当 海拔 在 3800~4800 m, 降 水 量 随 
海拔 的 增加 不 断 降低 ,在 海拔 高 于 4800 mm 后 , 除 
MLR 降水 量 随 着 海拔 上 升 而 小 幅 增加 外 ,其 余 降 水 
EHE FEBR, 
2.2 降 尺 度 模 型 评价 与 对 比 

为 了 对 降 尺 度 模型 进行 精度 验证 ,本 文 主要 利 
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图 3 研究 区 海拔 划分 (9) 及 降水 量 随 海拔 的 变化 (b) 


Fig. 3 The elevation division of the study area (a) and the change of precipitation with elevation (b) 


用 研究 区 域内 20 个 气象 站 点 的 年 均 . 季 均 .月 降水 
量 进行 对 比 验证 分 析 。 运 用 ArcGIS 中 的 提取 分 析 
工具 ,提取 各 站 点 降 尺 度 后 的 年 、 季 、 月 数据 ,并 计 
算 CC RMSE Bias 3 种 评价 指标 。 
2.2.1 年 降水 量 评 价 ” 从 图 4 中 可 以 看 出 ,TRMM 原 
始 数据 与 观测 数据 的 CC Al RMSE 值 为 0.83 和 
69.59;MLR .PCSR 和 Kriging 降 尺度 数据 与 观测 数据 
的 CC 值 分 别 为 0.86、0.85 和 0.85,RMSE 值 分 别 为 
73.39 .65.61 和 65.19。 从 这 个 结果 来 看 ,3 种 降 尺 度 
数据 得 到 的 降水 量 与 观测 数据 的 一 致 性 均 优 于 原 
台 TRMM 数据 ,其 中 MLR 数据 的 CC 值 最 优 。PCSR 
与 Kriging 数 据 的 RMSE 值 优 于 TRMM ,但 MLR 数据 
差 于 TRMM 数 据 。 此 外 ,Kriging 数据 的 Bias 值 优 于 
TRMM 原始 数据 ,PCSR 数据 没有 变化 ,而 MLR 则 差 
于 TRMM 原始 数据 。 综 合 来 看 ,在 年 尺度 上 基于 
Kriging 方 法 得 到 的 降水 数据 精度 最 优 ,其 次 为 PC- 
SR 数据 ,MLR 数据 精度 最 差 。 
为 了 讨论 TRMM 及 其 降 尺 度 结果 空间 分 布 与 
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气象 站 实际 降水 空间 分 布 的 一 致 性 ,以 气象 站 点 年 
均 降水 量 为 参考 ,计算 了 各 站 点 的 CC、RMSE 和 Bi- 
as, Ff LA 100.7°E 将 其 划分 为 西部 地 区 和 东部 地 区 ， 
以 此 得 到 年 尺度 上 东西 部 地 区 的 精度 评价 箱 线 图 
对 比 (图 5 和 图 6)。 在 空间 上 ,TRMM 及 3 种 降 尺 度 
数据 均 显 示 东 部 地 区 的 CC 值 较 高 ,RMSE 值 较 小 ， 
表明 该 地 区 格 点 数据 与 实测 降水 更 一 致 ,数据 的 误 
差 在 此 区 域 更 小 ;网 5 中 的 Bias 显示 在 研究 区 西 侧 
的 各 数据 均 低 估 了 降水 ,研究 区 东 侧 格 点 数据 高 估 
了 大 多 数 站 点 的 降水 。 由 图 6 结果 表明 ,TRMM 及 3 
种 降 尺 度数 据 的 CC 值 在 东部 地 区 均 高 于 西部 地 
区 ,RMSE 和 |Bias| 值 在 东部 地 区 低 于 西部 地 区 。 综 
上 所 述 ,TRMM 及 3 种 降 尺 度数 据 在 东部 地 区 精度 
最 优 。 

2.2.2 季 降 水 量 评价 ”由 图 7 可 知 , 在 数据 精度 方 
面 ,Kriging 数 据 在 春季 的 CC、RMSE 表现 最 好 ;从 夏 
季 来 看 , MLR 数据 的 CC 值 表现 最 好 ,但 RMSE K Bi- 
as 均 最 差 ,PCSR 与 Kriging 数 据 在 RMSE 与 Bias 值 表 
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图 4 实测 降水 量 与 TRMM 数 据 及 其 降 尺 度数 据 年 均 散 点 图 


Fig. 4 Annual average scatter plot of measured precipitation and TRMM data and its downscaled data 
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图 5 年 尺度 上 TRMM .MLR .PCSR 和 Kriging 的 CC(a~d)、PCSR(e~h) 和 Bias(i~]) 的 空间 分 布 
Fig. 5 Spatial distribution of CC (a-d), PCSR (e-h) and Bias (1-1) of TRMM, MLR, PCSR and Kriging on yearly scale 
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6 年 尺度 上 TRMM 及 3 种 降 尺 度数 据 在 东西 部 地 区 的 精度 评价 箱 线 图 对 比 


Fig.6 Comparison of boxplots for accuracy evaluation of TRMM and three types of downscaled data in the 


eastern and western regions on the annual scale 


现 均 较 好 ;秋季 的 RMSE 和 Bias 值 中 PCSR 数据 表现 
最 好 ,但 其 CC 值 较 差 于 Kriging 数 据 ;TRMM 数据 在 
冬季 的 CC、RMSE 和 Bias 表现 最 好 , Kriging 数据 表 
现 最 差 。 此 外 ,冬季 的 |Bias| 值 均 为 正 值 且 较 大 , 反 
映 出 4 种 数据 在 冬季 偏离 程度 均 较 大 ,日 均 大 于 实 
测 降水 量 。 以 往 研究 表明 ,在 站 点 尺度 上 ,普遍 认 


位 置 .三 散发 和 降水 类 型 等 ”“”。 其 中 ,降水 类 型 按 
照 气温 划分 ,主要 可 分 为 液态 降水 、 混 合 降水 和 固 
态 降 水 ,在 寒冷 的 青海 湖 及 周边 地 区 固态 降水 占有 
很 大 比重 ,而 雨量 简 对 固态 降水 的 捕捉 率 最 差 ” ， 
造成 该 地 区 冬季 的 实测 降水 量 低 于 实际 降水 量 。 
因此 ,在 冬季 固态 降水 观测 误差 本 身 比较 大 ,可 能 


为 雨量 简 获 取 的 降水 量 最 为 准确 ,然而 南 量 简 获 
取 的 降雨 量 受 很 多 因素 的 影响 ,如 自身 误差 ,安装 


是 冬季 TRMM 数 据 及 3 种 降 尺 度数 据 的 Bias 值 为 正 
值 且 比较 大 的 原因 。 综 合 来 看 ,春季 的 Kriging 数 据 
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图 7 TRMM 数 据 及 降 尺 度 季 均 数 据 精 度 对 比 
Fig.7 Accuracy comparison between TRMM data and downscaled quarterly average data 


精度 最 优 , 在 夏季 PCSR 与 Kriging 数 据 精 度 基 本 持 
平 ,在 秋季 PCSR 数据 最 好 ,而 冬季 的 TRMM 数据 表 
现 最 佳 。 

2.2.3 月 降水 量 评价 ”运用 提取 分 析 工 具 提 取 了 
2010—2019 年 各 个 站 点 12 个 月 的 TRMM 数据 、 
MLR, PCSR 及 Kriging 降 尺 度数 据 。 将 2010 一 2019 
年 的 月 气象 站 点 数据 降水 量 作为 自 变量 ,TRMM 数 
据 与 3 种 降 尺 度 降 水 数据 作为 因 变 量 ,得 到 散 点 图 
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(图 8)。 从 图 8 可 以 看 出 ,TRMM 数 据 与 实测 数据 相 
关系 数 CC 为 0.89, 均 方 根 误 差 RMSE 为 17.51, 相 对 
误差 Bias 为 -0.03;MLR 降 尺 度数 据 与 实测 数据 CC 
为 0.90, RMSE 为 17.36, Bias 为 -0.10;PCSR 降 尺度 
数据 与 实测 数据 CC 为 0.89, RMSE X 17.21, Bias 
为 -0.03; Kriging 降 尺度 数据 与 实测 数据 CC 为 0.72， 
RMSE 为 27.75,Bias 为 -0.02。 从 相关 系数 CC 来 看 ， 
MLR 降 尺 度数 据 同 实测 降水 数据 相关 性 最 好 ,其 次 
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图 8 TRMM 原 始 月 降水 量 、 降 尺度 月 降水 量 与 实测 月 降水 量 散 点 图 


Fig. 8 Scatter plot of TRMM original monthly precipitation, downscaled monthly precipitation and measured monthly precipitation 
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图 9 各 站 点 海拔 与 格 点 数据 精度 评价 指标 统计 


Fig.9 Statistics of the accuracy evaluation indicators of the altitude and grid data of each station 


为 PCSR 45 TRMM 数据 ,Kriging 数 据 同 实测 降水 数 
据 相 关 性 最 差 。 从 均 方 根 误差 RMSE 来 看 ,TRMM 
数据 和 3 种 降 尺 度 结果 与 实测 降水 值 的 接近 程度 表 
现 为 PCSR>MLR>TRMM>Kriging。 从 相对 误差 Bias 
来 看 ,|Biasl| 表 现 为 Kriging<PCSR=TRMM<MLR。 此 
外 ,通过 各 散 点 数据 与 y=x 线 进行 对 比 ,发 现 可 以 看 
出 ,PCSR 与 MLR 数据 集中 程度 相近 ,其 次 为 TRMM 
数据 ,Kriging 数据 离散 程度 最 高 ,精度 最 差 。 综 合 
来 看 ,在 月 尺度 上 PCSR 数据 精度 较 好 ,Kriging 数 据 
精度 最 差 。 
2.3 海拔 对 TRMM 及 其 降 尺 度数 据 的 影响 分 析 
以 往 研 究 表明 ,海拔 是 影响 山地 降水 的 重要 因 
素 , 随 着 海拔 的 升 高 ,降水 量 一般 呈 上 升 的 趋势 ,但 
遥感 数据 往往 无 法 准确 表达 降水 的 变化 ”。 为 此 
将 各 站 点 海拔 按照 由 低 到 高 的 顺序 与 其 对 应 的 格 
点 数据 精度 评价 指标 进行 统计 (图 9)。 结 果 表 明 ， 
CC 值 与 RMSE 值 随 着 海拔 的 升 高 没有 明显 的 变化 
趋势 ,是 波动 较 大 ,说明 海拔 对 研究 区 内 TRMM 及 
其 降 尺 度数 据 的 影响 不 大 ; 随 着 海拔 的 升 高 ,Bias {EL 
由 正 转 负 ,整体 呈 下 降 趋 势 ,表明 TRMM 数据 随 着 
海拔 的 升 高 逐渐 出 现 低估 降水 的 现象 ,这 与 卢 新 玉 
SEO) 、 李 琼 等 ”研究 结果 是 一 致 的 。 综 上 所 述 ,在 
1330~3800 m 海 拔 范 围 内 ,人 研究 区 降水 量 随 着 海拔 
的 升 高 而 逐渐 增加 ,因此 在 降水 量 较 多 的 山区 对 流 
性 降水 占 总 降水 比率 较 高 ,而 TRMM 数据 微波 降水 
率 反 演 时 对 对 流 性 降水 有 一 定 程度 低 佑 ““*” ,导致 
降水 量 较 高 的 山区 出 现 低估 降水 的 现象 。 


3 结论 


本 文 以 青海 湖 流 域 及 周边 区 域 作 为 研究 区 域 ， 


选取 经 纬度 .坡度 、 坡 向 .DEM、NDVI 作 为 影响 因 
子 , 对 2010 一 2019 年 间 的 TRMM 数据 进行 降 尺 度 处 
理 , 通 过 构建 多 元 线性 回归 模型 (MLR). 主 成 分 - 逐 
步 回 归 模 型 (PCSR) 和 普通 克 里 金 插 值 法 (Kriging) 
进行 降 尺 度 研究 ,并 运用 CC、RMSE Bias 作为 评价 
指标 进行 精度 验证 和 对 比分 析 , 得 到 以 下 结论 : 

(1) TRMM 及 3 种 降 尺 度数 据 与 气象 站 点 数据 
在 青海 湖 流 域 及 周边 区 域 得 到 降水 空间 分 布 具 有 
一 致 性 。 年 均 及 春 、 夏 、 秋 季 的 降水 量 均 表 现 为 北 
部 高 ,西部 及 西北 部 低 ,而 冬季 降水 量 表现 为 南部 
与 西北 部 高 ,中 部 低 。 此 外 ,研究 区 降水 量 随 着 海 
拔 的 增高 ,以 3800 m 为 界 整体 上 呈现 先 升 高 后 降低 
的 趋势 。 

(2) 在 精度 评价 方面 ,PCSR 降 尺 度数 据 在 青海 
湖 流域 及 周边 区 域 表 现 更 为 准确 。 从 年 尺度 上 来 
看 ,Kriging 精 度 表 现 最 好 ,在 空间 上 TRMM 及 3 种 降 
尺度 数据 在 东部 地 区 精度 最 优 ; 在 季 尺 度 上 ,数据 
精度 表现 为 PCSR>Kriging>TRMM>MLR ;在 月 尺度 
上 ,PCSR 数据 精度 最 好 。 

(3) 人 研究 区 内 TRMM 及 3 种 降 尺 度数 据 受 海拔 
影响 较 小 ,但 随 着 海拔 的 升 高 ,TRMM 及 其 降 尺 度数 
据 逐 渐 出 现 低估 降水 的 现象 ,其 可 能 原因 在 于 降水 
的 低估 与 微波 降水 率 反 演 时 对 对 流 性 降水 的 低估 
有 关 。 
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Abstract: Using multiple linear regression (MLR), principal component stepwise regression (PCSR), and 
Kriging, the TRMM 3B43 precipitation data in the Qinghai Lake Basin and surrounding areas with a resolution of 
0.25° were downscaled to a resolution of 0.01°. The measured precipitation data of 20 meteorological stations in 
the study area were selected, and the correlation coefficient, root mean square error, and relative deviation (Bias) 
were used to evaluate the downscaling results. Downscaling methods for the study area. The results show a 
consistent spatial distribution of precipitation in the study area based on TRMM and the three downscaling 
methods. The annual average precipitation and the three seasons of spring, summer, and autumn are all high in 
the north, low in the west and northwest, and winter precipitation. The performance was high in the south and 
northwest and low in the middle. With increased altitude, the precipitation in the study area showed an overall 
trend of first increasing and then decreasing with 3800 meters as the boundary. The results of the accuracy 
evaluation show that the Kriging accuracy on the annual scale has the best performance. On the spatial scale, the 
TRMM and the three downscale data have the best accuracy in the eastern region. On the quarterly scale, the data 
precision is PCSR > Kriging > TRMM > MLR. On the monthly scale, the PCSR data accuracy is the best. The 
effect of altitude on the TRMM and the three downscaling data in the study area is small. However, with 
increasing altitude, the remote sensing data gradually underestimates the precipitation phenomenon, which may 
be due to the underestimation of precipitation and convection during microwave precipitation rate inversion 
underestimation of precipitation. From the comprehensive precipitation spatial distribution consistency analysis 
and precision evaluation, it is considered that PCSR is the most suitable downscaling method for TRMM 3B43 
precipitation data in the Qinghai Lake Basin and surrounding areas. 

Keywords: TRMM 3B43; downscaling; multiple linear regression; principal component stepwise regression; 
kriging; Qinghai Lake Basin 


